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Abstract 

In this paper we extend previous hydrodynamic equations, governing the motion of Bose- 
Einstein-condensed fluids, to include temperature effects. This allows us to analyze some differ- 
ences between a normal fluid and a Bose-Einstein-condensed one. We show that, in close analogy 
with superfluid He-4, a Bose-Einstein-condensed fluid exhibits the mechanocaloric and thermome- 
chanical effects. In our approach we can explain both effects without using the hypothesis that 
the Bose-Einstein-condensed fluid has zero entropy. Such ideas could be investigated in existing 
experiments. 
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I. INTRODUCTION 



The achievement of Bose-Einstein condensation (BEC) in dilute atomic gases [1] repre- 
sents the estabhshment of several new exciting possibilities. The investigation of degenerate 
quantum gases provides, for instance, a good testing ground for innumerable many-body the- 
ories developed along the past few decades. Those theories were mainly developed within 
the context of explaining the striking properties of Helium-4 below the A-point. Among 
other features, this fluid behaves as having zero viscosity (a superfluid) and apparently also 
as having zero entropy [2]. 

The zero entropy hypothesis seemed to be the only explanation for two other (although 
related) remarkable properties of superfluid He-4. The flrst notable property is the called 
thermomechanical effect (or the "fountain" effect), where a temperature gradient across 
the fluid produces motion of matter. The force responsible for imparting the motion to 
the fluid is a "thermoforce" (a term used flrst by London [3]). On the other hand, the 
motion of the superfluid also produces temperature gradients. This second effect is the so 
called mechanocaloric effect. Both effects were very important to understand the behavior 
of superfluids, giving valuable hints about the nature of the processes involved. 

In earlier times the flrst proposal of London and Tisza [3, 4] , considering the superfluidity 
as a macroscopic manifestation of BEC, was somehow overruled by the good agreement 
of the hydrodynamic theory of Landau [5] to the flrst experiments measuring the second 
sound velocity [6]. Just later, after much of the development of the many-body theory, 
specially by Bogoliubov [7], it became clear that BEC was in fact behind the appearance of 
superfluidity. Even then, the hydrodynamic formalism has been quite useful for calculating 
and understanding the He-4 superfluid. Recently, after the experimental achievement of 
quantum degenerate gases, the formalisms developed initially to understand superfluidity of 
Helium have been extended to explore the dilute atomic gases as well [8]. A full historical 
account of these developments, although exciting, is much beyond the scope of this present 
paper. The interested reader, however, could look for instance at reference [9] and the 
references therein. 

In this paper we deal with some aspects of the hydrodynamics of Bose-Einstein condensed 
phase. In particular we show that both effects (the mechanocaloric and the thermomechan- 
ical effects) are always expected to happen in any Bose-Einstein condensate. 
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The way we deal with these two phenomena are derived from the same equation, showing 
the clear relation between them. Furthermore, we do not have to use the zero-entropy 
property for the superfiuid or the BEC system to explain both effects. 

The presentation of this paper is arranged as follows. In section 11 we introduce the hy- 
drodynamic quantization approach to the Bose-Einstein condensation of a system of charged 
spin-0 bosons under the action of external fields. A chemical potential is introduced later 
in order to discuss the equilibrium conditions for this system. In section III we present the 
equations of motion and define averages over the ensemble. 

In section IV we write the hydrodynamic equations of Bose-Einstein condensates under 
the action of external fields. These equations do not take into account temperature effects. 

The relevant equilibrium equation, a generalized Bernoulli equation, is introduced in 
section V. From this zero temperature equilibrium condition we then propose a natural 
extension that allows us to take into account also the temperature effects. 

The force equation and the explicit expression for the thermoforce is presented in sec- 
tion VI. Finally, the conclusions in section VII close the paper. 

II. HYDRODYNAMIC QUANTIZATION 

In earlier papers [10, 11] we have presented the hydrodynamic quantization approach 
and some aspects of the fluidity of Bose-Einstein condensed system. In this section we shall 
analyze these equations for the quite general and important case of a condensate in the 
presence of electric and magnetic fields and when the number of particles varies. In order 
to do so we introduce a chemical potential. Furthermore, the introduction of a chemical 
potential is also relevant in order to study the equilibrium conditions for condensation under 
the action of external fields [12]. 

We shall discuss here the field theoretic approach aimed at the description of a system 
of scalar Bose particles. Within the field theoretic approach, charged spin bosons are 
associated to a complex scalar field iIj{x). In what follows, x is a short hand notation for 
a time and space point {x = {t,x)). 

The quantization of the theory is carried out by imposing equal-time commutation re- 
lations between the field and its canonically conjugated variable followed by a typical 
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quantization scheme implemented with the introduction of the Fock space. This method is 
especially useful when we deal with systems whose number of particles is well defined, for 
example in the description of scattering processes. One can use, however, an alternative pro- 
cedure: the so-called hydrodynamic quantization of the field. In this scheme we use a new 
set of variables p{x) and ip{x), which are more convenient in the description of properties 
related to Bose-Einstein condensation. These new variables are defined as: 

p{x) — iIj*{x) ip{x) 

(II.l) 

ij{x) = ^/p{x) e^^(^) , 

denominated density and phase variables, respectively. They are canonically conjugated to 
each other, as we shall see below. For this reason we impose the following commutation 
relations among the density operator and the phase operator 

[p(x), if{x')] = i5(x-x') . (II.2) 

The quantization method, based on commutation relations among density and phase op- 
erators, is called a hydrodynamic quantization, and its application in superfiuidity, where 
Bose-Einstein condensation is known to occur, was proposed by Landau [12]. The commu- 
tation relation in equation (II. 2), in conjunction with an explicit representation for these 
operators, requires a departure from more usual quantization approaches. 

An explicit representation of the algebra of operators (II. 2) can be realized by making 
use of the so-called density representation. In this representation the density operator is a 
classical c-number. That is, 

p{x) = p{x) (II.3) 
and the phase operator is represented by the operator 




(11.4) 



Within the density representation the states vectors are represented as functionals of the 
density 

^ = m ■ (n.5) 

In Ref.[10] we give examples of how to construct wave functionals associated to the 
vacuum and Fock states. 
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The state functionals, as well as other relevant physical quantities, are written in terms 
of the density and phase operators (H.l), or in terms of integrals over space-time densities. 
We shall give below several examples that will be used through this paper. We start with 
the Lagrangian and the Hamiltonian. 

In terms of the field ■i/', the classical Lagrangian density for the nonrelativistic scalar 
particles, is written as 

where is the interaction Hamiltonian density. We define the Hamiltonian density 7i as 
composed of two terms: 

n{x) = Y^li±^^^ + n'[r{x)^{x)] = K + n' . (h.t) 

^ III 

The first term in (H.7), containing derivatives of the field, is the kinetic energy term (K) 
and the second is the interaction Hamiltonian Ti^ , which contains no field derivative terms. 

Placing equation (II. 1) into equation (II. 6) one can see that the classical Lagrangian 
density can be written under the general form: 

C{x) = -p(x)^-n[x,p,cp] . (II.8) 

Thus showing that phase and density are canonically conjugated variables. 
The quantum action is defined as 



S = (,1.9) 



If the system is under the action of external magnetic fields (B = 'V x A) the kinetic 
term can be written, in terms of the variables p{x) and (p{x), as: 

k{p,^,a)^ pix) L L + V J (1,10) 

where e is the electric charge of the spin bosons. 

Whereas, in the presence of an external potential U (x), the interaction Hamiltonian takes 
the general form: 

n\p) = p{x) U{x) + p{x) s{p, x) = p{x) U{x) + Ti}""' (H.ll) 



where s{p,x) in (11.11) is the per-particle interaction internal energy density. If we assume 
binary interactions among the particles, £{p, x) is given by 

£{p^x) = ^Jdx'V{x-x')p{x') . (11.12) 

Let us consider the expression, in terms of density and phase, of other physically relevant 
quantities in the study of quantum fluids. We start with the expression for the momen- 
tum. The classical definition of momentum is the integral over the space of the momentum 

operator density {V{x)) of the field: 



p = j dxV{x) = \ j dx(p*{x) V (f{x) = j d-'xp{x)Wip{x) . (11.13) 



From (11.13) it follows that the per-particle momentum P[x) is given by 

P(f) = V¥?(a;) = V ( -ip-^ ) = my(f ) . (11.14) 



5p{x) 

In quantum theory there is, however, another relevant momentum. This quantum mo- 

1 

mentum {Pq{x)) is defined as - of the gradient of the log of the density. That is 

Pq{x) = V (in VpM) ^mVq{x) . (11.15) 

The quantum origin of Pg{x) defined in (11.15) is explained in Ref. [14]. 

The analysis of a hard sphere gas provides a good example for the understanding of 
both types of momenta [11]. The classical momentum P{x) is associated to the quantized 

— * 

vortices, whereas the quantum momentum Pg is more relevant in the understanding the 
dependence of the fluid density upon the distance of the vortex center [11]. 

In terms of the velocities defined in (11.13) and (11.14) one can write, by using (11.10) and 
(11.11), the Hamiltonian density 7i as 

n[p, ^] = p{x) - \v\x) + V^{x)\ + p{x) [U{x) + e(p, x)] . (11.16) 

If the bosonic system is under the action of an external magnetic field the above expression 
remains valid for V{x) defined by 

V{x) = ""^"^ - ^^^^"^ . (11.17) 
m 
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The interaction of the particles with external fields is carried out, in field theory, as 
usual. For the coupling with external magnetic fields {B — V x A) we use the minimum 
substitution V ^ A — ieA. The potential U{x) takes into account the interaction with 
external electric fields as well as external gravitational fields. 

All the above definitions are relevant in the understanding of fiuidity aspects of Bose- 
Einstein condensed systems. 

In the following section we shall analyze the dynamical equations. 



III. EQUATIONS OF MOTION AND ENSEMBLE AVERAGES 

The time evolution of a physical quantity represented by the operator 0{x) is given by 
the Heisenberg equation of motion, 

^^i[H,0{x)] . (III.1) 

The Hamiltonian operator in (III.l) is the integral of the Hamiltonian density 



J dxn[(p,p] . (III.2) 

Within the density representation the equations of motion of the density and phase 
operators are: 

dip{x) 5H 



dt ^pipc) 

dpjx) _ SH 
dt 5(p{x) 



(III.3) 

(1II.4) 



These equations follow also from (II. 8) by taking p and (f as independent dynamical 
variables. 

By taking the gradient of (III. 3) we define an yet new equation that we name the force 
equation: 

^ = . (111.6) 

The force equation gives the rate of change of the per-particle momentum and it is the 
relevant equation in the understanding of superfiuidity and the thermoforce. As we shall 
see, the thermoforce is a new type of force that arises in a BEC system. 



The local force operator F{x) is, formally, written as 

One of the interesting features of the hydro dynamic quantization is that the quantum 
equation of motion resembles that of a classical fluid (from this fact derives the "hydrody- 
namic" name). Using Hamiltonian (11.16), the time evolution equations are 

dpjx) ^ _^ _ p{x) P{x) 
dt m ' 

f = {£} + - ^"(^j - ^ (£+ ^ ^ 

Where h{x) in (111.7) and (III. 9) is the per-particle enthalpy 

SW''^ P(t) 
h{x) ^ ^ =e(p,x) + 44. (III.IO) 
dp{x) p{x) 

The pressure P{x) in (III.IO) is, as can be inferred from (III.IO) and (11.11), given by 



p{x) 



J dx'p{x')^^e{p{x'),x') . (III.ll) 



— * 

The term Fe^t{x) in equation(III.9) is the force exerted on each of the charged particles 
of the system as a result of the external potential U{x), 

F,.t(x) = -VU(x) . (III.12) 

Equations (III.7)-(III.9) are vahd as time evolution equations of operators. In order to 
deal with classical c-numbers one can consider expectation values and averages over the 
ensemble. 

Any of the quantum versions of the equations of motion (III.7)-(III.9), are obviously 
valid when one considers expectation values. Defining the trace of an operator O as 

T,0{x) = J2{i^\0{x)\ij) , (III.13) 
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we can write the following equations: 



^(Tr^x)) = -tJ^] , (III.14) 



dt \5p{x 

i-(Trp{x)) = Tr(^) . (III.15) 



dt \S(p{x) 

Ehrcnfest theorems can be defined as averages over the ensemble of quantum equations. 
In order to define averages over the ensemble, we define the quantum-mechanical partition 
function of a Bose system. 

For the quantum action S, given by (II. 9), we define the partition function as 

Z = Tr e'^ . (III. 16) 

This is just an extension of the usual finite temperature (T) definition of the partition 
function: 

Z^Tr {e-^^ . (III.17) 
The average over the ensemble of a physical quantity 0{x) is usually defined by 

(0(x)) = ^Tr{e^'0{x)] . (III.18) 

Taking averages over the ensemble, one writes: 

d{ip{x)) I m 



dt \ 5p{x) I 



(111.19) 



m-)) - + / J« ) . (111.20) 



dt \ S(p{x) 

Averages over the ensemble are relevant in the context of statistical mechanics. In the 
next section we shall consider matrix elements of the equations of motion under the form 
(III.14) and (III.15). 



IV. HYDRODYNAMIC EQUATIONS FOR BOSE-EINSTEIN CONDENSED 
STATES 

One expects some basic distinctions, at the level of states, between a condensed system 
and a normal one. We have proposed, in Ref. [10], that the basic distinction can be traced 
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back to special properties of the wave functional associated to BEC states. We have proposed 
that the wave functional of BEC systems is endowed with two distinctive properties [10, 11]: 

Property (1) of the condensate wave functional. 

The wave function of the system associated to ■0c[p] [H, 15] is a product of wave functions 

TV 

tpcixi . . .Xn) =Y\_i'c{Xi) ■ (IV.l) 

i=l 

The wave function i^dx) in (ULl) is, by definition, the condensate wave function which 
will be written as 



M^) = y^e^^=(^) . (IV.2) 

The variables Pc{x) and fdx) will be identified as the density and phase of the conden- 
sate. 

Property (2) of the condensate wave functional. 

The condensate wave functional describes coherent states for which the following factor- 
ization property holds true: 



and 



M^)^ Vp;Re'^=(^) (IV.3) 



.^j^,^^ = Mx.) . ■ ■ M^n) . (IV.4) 

Property (1) is an obvious requirement in order that the wave functional be associated to 
a condensed state. Whereas property (2) allows us to identify the condensate wave function 
with the order parameter of the phase transition [16, 17]. 

The factorization property (111.4) is the analog of the Off Diagonal Long Range Order 
(ODLRO) introduced by Penrose and Onsager [18] in the context of superfluidity of He-4. 

In References [10] and [11] we have shown that the wave functional 

^^p] = e^^(^^'"'^=(^)'^^ (IV.5) 

exhibits the properties above mentioned properties (1) and (2). Thus providing an explicit 
example of such wave functionals. 
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The relevance of the factorization property in the understanding of superfluidity has been 
already emphasized by Anderson in Ref. [16]. By taking averages over states satisfying 
(IV.4), we conclude that the wave function of the condensate, defined in (IV. 3), satisfy the 
equations: 



dip^ix) ( m 



dt \5p{x) 
dpc{x) ( 5H 



dt \5(p{x 



ip=<Pc{x) 



Since a c-number field il){x) is the wave function in the x-representation of a state , 
equations (IV.6) and (IV. 7) specify the equation for the condensate wave function. They 
are c-number equations. 

By taking the functional derivatives of Ti. given by (II. 2), now expressed as a functional 
of p and we obtain the following set of equations: 

^ + V • Ux) = (IV.8) 
Where we have used for V, Vq and h{x) the definition (II.7), (11.15) and (III.IO). The 

— * — * 

density currents Jc and Jq are given by: 

Jq{x)^Pc{x) Vq{x) (IV.IO) 

J,{x)=p,{x)V,{x) . (IV.ll) 
Taking the gradient of (III. 9) we get the force equation: 

^ + m(v,-V^V, + m(Vq-V^Vq = eV,AB-V i^U + h+ ^V-pF,| . (IV.12) 

In the presence of an external magnetic field the motion is necessarily rotational since 
the rotational of the condensate velocity is related to the magnetic field by 

V xV^ = - — V X A = - — B . (IV.13) 
m m 
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In the absence of a magnetic field the motion of the condensate is irrotational. 

Equation (III. 8) is continuity equation which we have shown to be true for the condensed 
component of the system. 

As we shall see in the next section, equation (111.9) corresponds to a generalized Bernoulli 
equation. We refer to it as the equilibrium equation. It gives the condition of equilibrium 
for a Bose- Einstein condensed fluid under the action of external fields. 

Equation (III. 13) is the dynamic equation satisfied by the fluid in motion under external 
fields. We shall see that it is an extension of one of Anderson's equations in Ref. [16]. 

This allows us to conclude that the dynamical equations for the density and phase of the 
wave function of the condensate, are (IV. 8), (IV. 9) and (IV. 12). As a result of these equa- 
tions we can conclude, on quite general grounds, that the condensed fluid satisfy continuity 
equation and that in the absence of magnetic fields the motion of the fluid is irrotational. 



V. EQUILIBRIUM EQUATION 

We show in this section that, for stationary states, equation (IV.9) gives the equilibrium 
condition when the system is under the action of external fields and that this equilibrium 
condition is a generalized Bernoulli equation. 

As usual, we assume that the equilibrium condition for a system under the action of 
external fields is the chemical potential (fj) to be constant: 

= /xo • (V.l) 

Remembering that the field tp is the wave-function of a state i/j in the f-representation, 
the field associated to a stationary state is 

iP{x,t) = e<'^"*+'^^(^")) (V.2) 

where fiQ is a constant. 

For a stationary state the equilibrium equation, for the phase given by (V.2), is 

^0 = h{x) + U{x) + ^ + - Vg\x) - -^^j^ ■ (V.3) 

Taking now into account the explicit expression for h{x), we write 



12 



, , ,x P{x) m -2 X TT, \ ■rr2i \ V • Jn(x) .^^ ^. 

lio = s{x,pc{x)) + + - V,\x) + U{x) + - V^\x) - ■ (V.4) 

In order to show that /iq is the chemical potential we multiply equation (V.4) by Pc{x) 
and integrate over x . The result can be written as 

fXoNc^ E, + j P{x)d^x (V.5) 

where Nc in (V.5) is the number of particles in the condensate 

N, = j d^xp,{x) (V.6) 

and Ec is the condensate energy defined as the integral of the density Hamiltonian: 

E, = J d'xn[pc,V^c] ■ (V.7) 

Expression (V.5) shows that Pq is the chemical potential since PqNc is the total enthalpy 
of the system [12] . 

For a uniform condensate, the Bogoliubov condensate [20], 

Pc{x) = pQc ■ (V.8) 

The quantum velocity vanishes so that for a uniform condensate, equation (IV. 3) is 
Bernoulli's equation: 

po = e{x, poe{x)) + m ^ + ^ ^2(0;) + U{x) . (V.9) 

POc ^ 

Since the first four terms of equation (IV. 4) corresponds to the usual Bernoulli equation, 
one can predict that deviations from the usual fiuid is expected as a result of the quantum 
velocity term. This prediction deserves to be tested experimentally. 

We notice that for a non self interacting system 

h{x) = . (V.IO) 

Such non self interacting fluid, the fluid described by (IV. 2) and (IV. 3), is a Mahdclung fluid 
[19]. Therefore, neglecting self interaction of particles, the BEC is a Mahdelung fluid. 

One can extend some of the previous equations to finite temperatures. In order to do so 
we just recall that the equilibrium condition for a fiuid under external fields takes the form 
[12]: 

p^g{P,T)+p' (V.ll) 
13 



where g{P, T) is the per-particle Gibbs energy in the absence of the external fields and fj,' 
takes into account the external fields. 

Remembering also that at zero temperatures the per-particle Gibbs free energy is equal 
to the per-particle enthalpy 

g{P, 0) = e {x, poc) + m = h{x) (V.12) 

we can write equation (IV. 6) under the following form 

/.o = g{P, 0) + - V^{x) + U{x) . (V.13) 

The natural extension of (V.5) to finite temperatures, in view of (V.13), is 

/^o=^(^,r) + ^^^^ + c/(x) 

=e ix, po) + m -Ts(T) + m^ + U(x) . (V. 14) 

Po 2 

Notice that for non stationary states the generalization of (IV. 11) to finite temperatures 
should be, from (V.2) and (V.14) 

^ = ^(P,r)+m^+C/(x). (V.15) 

The interesting aspect of (V.14) is that one can imagine an Entropy Filter. That is, 
one can devise experiments by means of which a low entropy component is separated from 
a large entropy component of the BEG fiuid. This would be similar to what happens in 
superfluid IIe-4. 

In fact, one can think of two experimental set ups in order to check the validity of the 
predictions made here. These experiments are analogous to the ones done with superfiuid 
He-4. Let us consider two containers (^4 and B) of a Bose-Einstein condensed fiuid in which 
we keep the temperature initially constant on each side and keep the density constant. At 
this point we consider no external fields. Suppose that these containers are now physically 
connected by a tube, in which we can have a flow of the condensate. Since 

dg ^ -^dP-sdT (V.16) 
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it follows from (V.14), that 



-dP-sdT ^ - d 
P 



2 



) 



(V.17) 




The temperature difference giving rise to a pressure difference is the thermomechanical 
effect. Equation (V.18) is London relation [3]. 

If the containers discussed above are kept at constant pressure, then from (V.16) it follows 
that 



As a consequence of (V.19) one can see that if there is mass flow in one of the containers, 
it will be colder than the other one. That is, the mechanocaloric effect is expected to happen 
in any BEC system. 

We can see, from (V.19), that in any BEC system a gradient temperature leads to motion 
of the fluid and vice- versa. It follows, also from (V.19) that, along a stream hne, the difference 
between the kinetic energy of the fluid at two points B and A, in each container, is: 



Prom equation (V.20) it follows that if the fluid flows from A to B {V{B) > V{A)), the 
container B will be warmer. On the other hand, if one of the containers is now considered 
as a thick tube the pressure dP will cause a fountain of Bose-Einstein condensate. We 
predict, in this way, the fountain effect for any BEC system. 

In superfluid He-4 these two phenomena are explained by assuming that the superfluid 
component is a zero-entropy fluid. We have seen here that this is not necessary. 

VI. FORCE EQUATION 




(V.19) 



mV^{B) _ mV^{A) 
2 2 



s{T) {Tb - Ta) 



(V.20) 



In this section we will show that the phenomenon of superfluidity for a BEC system is inti- 
mately connected to the two effects present in He-4, the thermomechanic and mechanocaloric 
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effects. Usually, they are considered as independent properties of the superfluid Hehum. We 
shall give, also in this section, formal expressions for the thermoforce. 
First of all we recall that 

dg = -^dP-s{T)dT . (VI. 1) 

Now we take the gradient of (V.14), and the dynamic equation for a nearly uniform fluid 
becomes 

— * 

dPc ( ^\ ^-t ^ ra 



^ +m(^K-Vjye = ^ext - -VP + eKxS + s(T)Vr. (VI.2) 

All terms in equation (VI.2) have a simple interpretation. The term in the left hand side 
of equation (VI.2) is the total derivative of P^.^ that is: 

DP, ( d ^ ^\ ^ , , 

^ U7 + ^c-V Pc. (VI.3) 



Dt \dt 

The first term in the right hand side of (VI.2) is the external force 

Pext = -VC/ . (VI.4) 

The second term in the right hand side of (VI.2) is the usual force that induces motion 
in the fluid as a result of pressure gradients 

Fi ^ - -VP . (VI.5) 
P 

The third term in the r.h.s. of (VI.2) is the Lorentz force due to the external magnetic 
field 

-^Lorentz — 

eVcXB. (VI.6) 

The last term in (VI.2) is our main contribution to the BEC motion. We propose that 
this term is the thermoforce: 

Pther = s(r)Vr. (VI.7) 

The presence of the thermoforce in equation (VI.2) means that there will be motion of the 
fluid as a result of temperatures gradients. The motion is towards the greatest temperatures. 
Putting all this together, we write 

^ = Fext + ^her + ^VP + eKxP . (VI.8) 
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Except for the thermoforce, equation (VI. 8) is just the Euler equation for an ideal fluid. 
Under similar assumptions made here, Anderson derived the usual Euler equation for He-4 
[16]. 

The conclusion, therefore, is that the same phenomenon (Bose-Einstein condensation) is 
responsible for all the striking properties of superfiuid He-4: the superfluidity and the both 
related thermomechanical and mechanocaloric effects. In fact, in principle, we expect that 
any BEC system should present those phenomena. At least one of them, has already been 
experimentally confirmed [21]. 

VII. CONCLUSIONS 

We have shown in this paper that the thermomechanical and mechanocaloric effects are 
always expected to occur in a Bose-Einstein condensate. Besides the understanding of these 
effects, we have shown that they can be predicted by using the same equation that leads to 
Euler equation (the equation that wc have named force equation). 

We have proposed a generalized Bernoulli equation that would lead, for non-uniform 
density fluids, to a departure from the usual Bernoulli equation. The generalized Bernoulli 
equation is valid for stationary BEC states. For time-independent states the chemical po- 
tential is zero, in this case the Bose-Einstein condensed state is characterized as a zero Gibbs 
energy fluid. 

Another important property of a superfiuid is the fact that it is irrotational [22]. As a 
consequence of the presented calculation, we have deduced that this is not generally the 
case. It only happens in the absence of external magnetic fields. 

We have derived an expression for the thermoforce. It is interesting to note that the 
same expression was also found by London, [3], using a different approach. Experimental 
evidence for the specific form (VI. 7) is, to our knowledge, missing. However, it predicts that 
as the temperature is lowered the thermoforce tends to zero (since lim s(T) — > 0). This 
fact is quite well known in superfiuid Helium-4. The thermoforce becomes negligible for 
temperatures below 0.6°K [3]. 

Another important contribution of this present paper is that in order to explain some 
aspects of the BEC phenomenon, there is no need of the zero entropy fiuid hypothesis. 

Finally we would like to mention that some of our specific quantitative predictions, con- 
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tained in the presented equations can be tested using present experimental conditions. 
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